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A scientific model need not be a passive and static descriptor of its subject. If the subject is 
aflfected by the model, the model must be updated to explain its affected subject. In this study, 
two models regarding the dynamics of model aware systems are presented. The first explores the 
behavior of "prediction seeking" (PSP) and "prediction avoiding" (PSP) populations under the 
influence of a model that describes them. The second explores the publishing behavior of a group 
of experimentalists coupled to a model by means of confirmation bias. It is found that model aware 
systems can exhibit convergent random or oscillatory behavior and display universal 1/f noise. 
A numerical simulation of the physical experimentalists is compared with actual publications of 
neutron life time and A mass measurements and is in good quantitative agreement. 



INTRODUCTION 

The notion that the act of observation can alter the 
observed system is familiar to any contemporary scien- 
tist. Less familiar is the notion that a model can al- 
ter the system it aims to describe. Such model-system 
couplings may be substantial for cognitive [HIS], social[3l- 
[7], or economic [5HTU] systems in which the constituent 
agents have full or partial access to the model that rep- 
resents them. It is conceivable for example, that an anal- 
ysis of word frequencies in spoken or written language [11 j 
can alter the behavior of those who produce these words; 
an analysis of stock market [12] can alter the behavior 
of stock traders; and an analysis of fashion[T3j can alter 
the behavior of those who influence fashion. The "phys- 
ical sciences" do not seem to be free of model-subject 
interaction either, since experimenters are influenced by 
theory through confirmation bias[15] and theorists build 
new models based on these experiments. 

The interaction between subjects and their models has 
long been recognized in separate disciplines, and referred 
as "the looping effect" in philosophy, "the enlighten- 
ment effect" [3 in psychology, and "performativity" [8j or 
"virtualism" [5] in economics. The model aware behav- 
ior in the physical sciences has been recognized too, and 
rightfully renounced as "confirmation bias" '[14l[T5]. Un- 
fortunately, the cited considerations are very qualitative, 
and discipline-specific. Presently there exists no quan- 
titative and generally applicable theory of model aware- 
ness. Without such a theory, it is not possible to ask 
whether the description of a model aware system will 
converge to a fixed point, or perpetually change while 
manipulating its subject. Neither is it possible to ask 
how the behavior of a system changes with increasing 
model awareness, or what universal properties, if any, do 
model aware systems have. 

The present study aims to answer these questions, at 
least partially, by introducing two quantitative models of 
model awareness, cast as much as possible in a discipline- 
independent language. The first describes a population of 



prediction seeking or prediction avoiding agents who up- 
date their behavior at every time step depending on the 
current model that describes them. The second describes 
the behavior of a population of experimental physical sci- 
entists, who decide whether to publish or not depending 
on the proximity of their data to the current "model" . 
In both the physical and social case, populations pro- 
vide feedback to the models too, since models must be 
updated to explain the current state of the population. 
The outcomes are then compared against two behav- 
ioral experiments [H [5] qualitatively, and particle physics 
data|16| quantitatively. 

Even though real-life models usually involve elaborate 
verbal descriptions, and/or sophisticated mathematical 
machinery, in this study we will consider much simpler 
ones such as propositions regarding majority behavior, 
or averages of multiple publications regarding a physical 
quantity. As simple as they may be, it seems appropriate 
to refer to them as "models" since they are falsifiable 
descriptions with predictive power. 

MODELS 
Model Aware Social Populations 

There has been a number of agent-based herding/anti- 
herding approaches in the literature used to explain a 
wide variety of social and economical phenomena (the 
interested reader is referred to [H]). The present ap- 
proach is reminiscent of a Polya urn process [18], in the 
sense that we will be modifying an ensemble according 
to its sampled outcomes. 

We consider a population of ^ 1, in which an indi- 
vidual can be in either one of the states A or B. These 
states can represent any opinion, property and behavior. 
At any given time t the state of the population can be 
characterized by the fraction of < (/)(i) = At/N < 1 of 
individuals in state A. After each time step, a scientist 
performs a measurement by choosing n random Individ- 



2 



uals (1 ^ n ^ iV) out of TV and publishes whether the 
majority of n is A or i?. In connection with our intro- 
ductory remarks we take any statement regarding major- 
ity behavior as a model; one example may be "investors 
avoid risk" . 

When the model is published, each individual among 

N may become "aware" of the model with probabihty 
<C 1, and subsequently update his or her state. The 
state of those unaware of the model is assumed to remain 
unchanged. 

Two types of populations will be considered separately, 
defined in terms of how the aware individuals update 
their state: A prediction seeking population (PSP) is one 
in which the aware individuals align their state with the 
prediction of the published model. For example, if the 
publication reports that "most people are A", then the 
aware B's become A in the next step (the aware A's 
remain unchanged). A prediction avoiding population 
(PAP) is one in which individuals anti-align their state. 
For example, if the publication reports that "most peo- 
ple are A" , then the aware A's flip their state to B (the 
aware B's remain imchanged). 

Once the aware subpopulation updates its state, a sub- 
sequent scientist will be sampling, measuring and model- 
ing a population of different nature. We will consider the 
Markovian dynamics of (j){t) over many such iterations. 
Unlike a Polya process, our ensemble is fixed in size, the 
number of modified agents is nondeterministic, and the 
modification is done according to a majority rule. 

Let us first focus on a PSP. Given the state (j>{t) — 
Ai/N , the probability of having k people who are A's in 
a sample of n is given by the hypergeometric distribution, 

'At\ fN - At^ 

PA{k) = (1) 
n 

Without loss of generality, suppose n is odd. The prob- 
ability that the majority of the scientist's measurement 
sample is A is 

n 

P{MA\At) = l-P{MB\At)= J2 PA{k) (2) 

fc>n/2 

After the publication, the probability P{r) that r people 
will become aware of the result is 



P(r) 



/(l-g)^--- 



(3) 



Given that there exists r aware individuals, and given the 
majority of the sample is A, the probability that u > 
of them {u < r) to change into B is 

At\m- At^ 

,r — ul \ u 



P{u\r,MA) 



N 



(4) 



Thus, the probability that the population changes from 
At to Atj^i =At + u\s given by the transition matrix 



N 



Ta,a+u = J2 ^("1^' MA)PiMA\At)P{r) 

r>u 

^P{u\qN,MA)P{MA\At) 



(5) 



Similarly, the probability that u > individuals change 
in the other direction is 



N 



Taa-u = J2Piu\r,MB)P{MB\At)P{r) 

^ P{u\qN,MB)P{MB\At) 



(6) 



where this time, 

P(u|r,Ms) = 



- At 
u J \ r — u 

r 



(7) 



We immediately see that if At = N, P{MA\At) = 1, 
P{u\r, Ma) ~ 1 and Fatjv = 1- The same holds true if 
At = 0, and Fo,o = 1- Thus, A = 1 and A = are the 
absorbing states of the corresponding Markov Chain. 

In the case of PAP, regardless of what the majority in 
the selected sample is, there always is a nonzero prob- 
ability that the aware population will include some of 
individuals from the majority state, who in the next 
time step will flip. Any accidental trend due to finite 
sampling n is counter balanced in the next time step by 
the prediction avoiders, who consequently cause an anti- 
trend. Thus, a PAP state will fluctuate on average by g0 
(though {\imt-i. oo (l>{t)) exists; cf. below). The maximum 
expected change {umax) = Nq will occur when At = N 
or At = 0. 

Let us now consider the expectation value (</>). For a 

PSP the expectation value of the change (Ac/)) is q{l — 
(t){t)) UMa, and -q4>{i) li Mb- Therefore, 

(A</)) = g(l - cj,)P{MA\At) - qct>{t)P{MB\At) 

= q[P{MA\At) - cl>]. (8) 

Similarly, for a PAP (A(/)) is -q(l){t) if Ma and g(l-(/)(t)) 
a Mb. Therefore, 



(A^) = q[l - P{MA\At) - <!>]. 



(9) 



The function P{Ma\(I)) < (j) for cj) < 1/2, P{MA\(t)) = (j) 
for (j) = 1/2 and P{Ma\(I)) > </> for > 1/2. Thus, the 
(;&c = 1/2 is a first order equilibrium state for both PSP 
and PAP. For a PSP, {Ac/}) is positive for > 1/2 and 
(lim(_>oo 0(i)) = 1 if (/'(O) > 1/2. (A0) is negative for (p < 
1/2, and {liuit^oo <f>{t)) = if (/>(0) < 1/2. In contrast, 
for a PAP, (A(/)) is negative for (/> > 1/2 and positive 
for (j) < 1/2; hence {limt^oo 4'{t)) = 1/2 regardless of 
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FIG. 1: Fraction of individuals (p{t) having a 
certain quality A as a function of time t. 

Prediction seeking (blue) and prediction avoiding 
populations (red) evolve under the influence of an 
unbiased model (left) and a biassed e=0.45 model 
(right). Here, iV = 1 x 10*, n = 5 x 10^, q = 0.05 



FIG. 2: < A0(^) > as a function of 0. Prediction 
seeking (left) and prediction avoiding populations 
(right) in the presence of error e = (solid) and e = 0.45 
(dashed), e shifts the equilibrium point for both PSP 
and PAP, and causes regular oscillations in PAP. Plot 
parameters are iV = 1 x 10*, n = 5 x 10^, q ~ 0.05. 



the initial state. Note that on average, the variance 
within the population is monotonically decreasing for a 
PSP, and monotonically increasing for a PAP (Fig(l)). 

It is also interesting to study the effects of incorrect 
models. This could happen for example due to an asym- 
metry in identifying one of the states, or misinterpret- 
ing correctly measured data during "modeling" . Suppose 
that the scientist publishes (1 — e)(/) instead of cj), the ef- 
fect of which can be taken into account by modifying 
P{MA\At); 



P{MA\AtY = I - PiMelAt)' 

n 



(10) 



/c>n/[2(l-e)] 

By observing vs (j) (cf. eqn (|8]),([9| and Fig(2)), we 
see that e shifts the equilibrium point (f)c for both PAP 
and PSP. As a result, the PSP may now cross (f> = 0.5, 
and unlike the e = case the variation within a PSP pop- 
ulation need not monotonically decrease. The absorbing 
states of a PSP does not change, whereas for a PAP e 
simply shifts the value of (limt_j.oo 0) ■ 



Numerical simulations are carried out for N = 10*, 
n = 500, q — 0.05, reveal interesting features regarding 
fluctuations (Fig(l), Fig(3)): First, in a PAP, as the error 
e is gradually increased within < e < 0.5, the period 
of the small fluctuations in equilibrium dramatically in- 
crease and regularize (Fig(l)). The reason is evident in 
Fig(2); for e = 0, we have |A0((/., - S)\ = |A(0, + 5)|, 
and therefore the period of trend/anti-trend fluctua- 
tions is equal to one time step. On the other hand for 
< e < 0.5, the difference in A0(0c - S) and A((/)c + S) 
leads to slower rises following rapid drops. As we con- 
tinue to increase 0.5 < e < 1, eqnQ now has a stable 
point, and the PAP starts to display convergent behav- 
ior much like the PSP. Thus, one could say that an un- 
changing model of a PAP is only possible if the model is 
practically wrong. 

A more interesting feature is that despite their oppo- 
site nature, PSP's and PAP's both appear to exhibit a 
1// noise in their power spectrum P{f ) if g <C 1 (Fig(3)). 
The 1// noise is the signature of self-organized critical 
systems HI] and is ubiquitous in nature (see for eg. |2IH- 
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\T2\). We predict that the same should appear in mea- 
surements of social model aware populations. 




FIG. 3: The power spectrum of prediction 

avoiding (red) and prediction seeking (blue) 
populations. It can be seen that both populations 



obey l/u} statistics (green). Here N 
q = 10-3. 
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The Model Aware Behavior of Physical Scientists 

Ideally, separate measurements of a physical quantity 
typically fluctuates around a fixed value according to a 
normal distribution due to random measurement errors. 
However, since experimenters tend to be influenced by 
past measurements and established theoretical models, 
data points are highly correlated [15] and have visible sys- 
tematic trends over timejlB]. This may be because, if the 
outcome of an experiment is significantly different from 
an accepted model or past outcomes, the experimenters 
may be "improving" their setup by eliminating some of 
the systematic errors which otherwise average out to zero, 
or even by simply repeating the experiment. In turn, fu- 
ture models are built on experimental data that have 
already been influenced by past models. Thus, it seems 
that the presence of model aware behavior extends be- 
yond social sciences. 

We describe the behavior of a population of model 
aware physical scientists with the following simplified 
characteristics: An experimenter measures a physical 
quantity fj, at time t and obtains (j)(t) in order to publish 
the n*'' paper P„ on the subject. Because of the precision 
limitations of the apparatus, (f>(t) is normally distributed 
around /i with standard deviation cq. After taking his 
measurement, the experimenter compares his result with 
the average Rq{n) = q~^Y^l=n-q+i of q preceeding 
published measurements, with standard deviation Sq{n). 

In this case we define the "model" to be Rq{n) (or 
any set of verbal or mathematical axioms that reproduce 



it). If (j){t) ± eo does not overlap with Rq{n) ± Sq{n), 
the experimentalist suspects there might be a mistake in 
the setup and decides to repeat the experiment at t + 
1. If on the other hand, if {P)q ± Sq overlaps with the 
"model", he publishes P„+i = (/)(t). It is assumed that /Lt 
is independent of time, and that non-random errors are 
negligible. 

The random variable Pn+i depends on the past q pub- 
lications {Pni Pn-i, ■ ■ ■ 1 Pn-q+i} ^nd is distributed ac- 
cording to a "partial Gaussian" function. 



/•Fl.q{n)+Sq+eo 

(Pn+i) = A{e,fi,s,R) pf^,eo{p)dp (11) 

J Rq(n) — Sq—ea 

where f^,eo(p) is the usual Gaussian distribution with 
mean and standard deviation cq. Since the possible 
publications are limited by the scientist's confirmation 
bias, the normalization constant A{e, ^, s, R) is defined 

We can immediately see that when Sq ^ eg, the Ihs 

/I regardless of P„; 



of (11) is precisely fi and (Pn+i) 



in other words, a lack of scientific consensus speeds up 
model convergence. The integral and A can be calculated 
exactly. 



-{eo+f^-R+sf/(2el) 



(eo-p+-R+s)V(2eo) 




(P. 



Note that if R{n) = fi we get (Pn+i) = M; ^^^d then 
R{n + 1) = R{n). Also, note that if R{n) < fi we have 
(P„+i) > P(n), and as a result R{n -I- 1) > R{n). On 
the other hand if R{n) > /i we have (Pn+i) < R{n) and 
as a result R{n -I- 1) < R{n). Thus P = ^ is a stable 
equilibrium point. Furthermore, the quantity (Pn+i) — 
R.fi is linear in — (P(n) — fi) up to third order in R{n) — /x, 
hence for values R{n) ^ Pn ~ fi we have exponential 
convergence (cf. Fig(4)). 

Fig(4) shows the outcomes of exact numerical simu- 
lations and compares them with historical publications. 
The first two data points are taken from actual particle 
physics experiments [TB] as an input, and the rest of the 
data points are iterated. The short time behavior is char- 
acterized by tight clusters separated by abrupt jumps. 
The long time behavior approximately fits an exponen- 
tial relaxation P„(t) = fj, + [fi — eo(0)]e-"/'^ curve. The 
convergence time T{q) determined from a least square 
fit, is a random variable with mean and standard devia- 
tion plotted as a function of q (Fig(5)). The convergence 
times increase and diversify with increasing model aware- 
ness of the experimenters. We conclude that model aware 
experimenters considerably slow down scientific progress. 
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FIG. 4: Published values of a physical quantity P„ at publication number n. Here, the Average values of 
publications of a simulated population of physical scientists (circles) is compared with actual |16j particle physics 

measurements (boxes) of neutron life-time (left) and A mass (right) normalized by their last value, g = 4 is chosen 
for the neutron measurements and q = 8 for the A measurements. The two initial points are taken from data and 

used as initial conditions in the simulation. 



FIG. 5: Average convergence time r as a function 
of model awareness q. Here two identical 
measurements with 15% error is taken as initial 
conditions. The error bars indicate the standard 
deviation in t. 



DISCUSSION 

To further connect model awareness to the experimen- 
tal literature, two recent behavioral studies [U |2j will be 
discussed within the PSP / PAP framework and will be 
compared with the simulation results. Due to lack of time 
dependent behavioral data, this comparison will have to 
be qualitative. Fortunately it will be possible to compare 
the outcomes of the simulations of the physical scientists 
quantitatively, since there exists records of historical par- 
ticle physics publications regarding single quantities |16j . 

The first behavioral study is one on the behavioral ef- 
fects of belief of testosterone administration and actual 



testosterone administration 2^. It was found that the 
women who believed that they were given the hormone 
behaved significantly more unfair than those who believed 
they were not given the hormone (regardless of whether 
or not they actually were) . Curiously the effect of actual 
testosterone was to increase fairness. 

Measuring the effect of belief of testosterone admin- 
istration relates to measuring the effect of the model of 
testosterone on its model aware subject. The model of 
testosterone is the (as it turns out, incorrect) statement 
that "testosterone causes people to be unfair" . The sub- 
jects of this study were familiar with this model, and 
behaved in accordance with its prediction. However, the 
study not only measures the effect of the testosterone 
model, it also alters it. Therefore if one repeats the same 
experiment in a few years, one may find that the effect of 
model to be the exact opposite, since if q is large enough 
some participants may be familiar to the altered model 
of testosterone. 

The second behavioral study is one on the effect of so- 
phistication (defined as the knowledge of psychology and 
statistics) on the ability of generating random numbers 
[T]. Here it was found that sophistication can dramati- 
cally reduce -if not entirely eliminate- the many nonran- 
dom trends unsophisticated people tend to display. Some 
such trends include avoiding repetitive triplets (such as 
8,8,8) or favoring descending sequences (such as 5,4,3) 
over ascending ones (such as 3,4,5). Here too, the author 
measures the influence of a model, which in this case, 
describes how people generate random numbers. Some 
of the people, (labeled "sophisticated") are aware of this 
model, and actively avoid its predictions. However, now 
that some people are aware of how people generate ran- 
dom numbers, the model must be updated to included 
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the behavior of this sophisticated fraction. Further- 
more, if q is large enough, upon repeating the experiment 
one might see the influence of "super-sophisticated" sub- 
jects, who avoid the trends that the mixed sophisticated- 
unsophisticated population display. 

Unfortunately both of these experiments have only two 
data points that are useful for us. The belief of testos- 
terone administration activates the present stage of a 
long model aware evolution in a PSP. Similarly the "so- 
phistication" variable in the random number study is a 
marker of a similar final state of a PAP. Thus one may 
only compare 0o to (j)t and Ao to At, where t is the 
present time. The short term behavior of the simula- 
tions qualitatively agree with both experiments [HE], the 
common features of which are (i) the model aware sub- 
jects are significantly different than the non model aware 
subjects and (ii) the non model aware subjects are more 
similar among themselves than model aware subjects. 

It is appropriate to represent the model aware subjects 
of the testosterone study (i.e. those who believe they re- 
ceived testosterone) by a PSP evolving under a strongly 
biased model. The model is taken to be biased because 
the statement that "testosterone makes people unfair" 
was demonstrated to be incorrect in the study. The sim- 
ulation agrees with feature (i) and for short times, fea- 
ture (ii). For long times, the simulated PSP eventually 
becomes less diverse than its starting state. Perhaps con- 
tinued publications of studies similar to [2] over the years, 
will reveal this time dependence. It is appropriate to rep- 
resent the model aware subjects of the random number 
study [T with a PAP evolving under an unbiased model. 
Here, too the simulation agrees with feature (i) and fea- 
ture (ii) (remember that (f> — 0.5 represents a maximally 
diverse population). 

The actual neutron half-life and A-mass measurements 
from 1960 to 2010 as reported in [16] are compared with 
the outcomes of the simulations, and is good agreement 
(Fig(4)). It is very interesting that other historical parti- 
cle physics data reported in [TB] fits reasonably well to the 
same exponential form. The average convergence time r 
depends on q. For example, since r is larger for A mass 
measurements compared to neutron life-time measure- 
ments, one can conclude that the scientists measuring A 
mass were more model aware. 



CONCLUSION 

It was shown that for physical scientists and prediction 
seeking populations, models and their subjects can co- 
evolve to a consistent state. Loosely speaking, for these 
systems model awareness sets a lower time limit to reach 
"truth" . For prediction avoiding populations, such a con- 
sistent state is possible only if the systematic error in the 
model is large enough (e > 0.5) to make the model practi- 
cally incorrect. For these systems, model awareness pre- 



vents convergence all together. This suggests that there 
may exist deterministic mechanical systems that are un- 
predictable due to their coupling to the predictor. The 
1 /f spectrum observed in populations of opposite nature 
suggest that there may be other universalities common 
to all model aware systems. 

While this is just a preliminary study, the simulations 
presented here demonstrate that the models and their 
subjects can be highly coupled and radically alter each 
other. Since the interaction of subjects with models seem 
to be present in a very diverse range of fields, the frame- 
work proposed was intentionally kept simple. This way, 
the theorists of different fields can add relevant system- 
specific details, and study the variants of the proposed 
model. Such variants could include interaction between 
individuals, noise, or variable model errors. 

Our study opens a wide range of additional questions 
that can be explored theoretically and experimentally: 
What is the influence of a time dependent model? Is it 
possible to construct a model that takes into account its 
own influence? For example, can a model predict its own 
acceptance in a community? Can the above considera- 
tions be generalized to more realistic models involving 
not mere numbers, but simulations or equations? How 
about model aware systems such as the stock market, 
where the subject is under the influence of multiple mod- 
els? Hopefully our work will inspire the scientiflc com- 
munity for a deeper exploration of model aware systems. 
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Appendix: Mathematica Codes 
(*The PSP/PAP Simulation*) 

Clear [u, p, n, q, i, average]; 

u = 10000; (*total population*) 

n = 500; (*experimenters sample population*) 

steps = 150;(*total time steps*) 

q = 0.01; (*awareness probability*) 

p = Table [0, {i, 1, steps}]; 

average = Table [0, {i, 1, steps - 1>] ; 

varicince = Table [0, {i, 1, steps}]; 

ensemble = Table [0, {i, 1, steps}]; 



Flip [number_] := Mod [number + 1, 2]; 
Flipper [p_, average_, q_] : = 

Table [RandomChoice[{l - q, q} -> {p[[i]]. 
If [Abs [average - p[[i]]] <= 0.5, 

Flip[p[[i]]], p[[i]]]}], {i, 1, u}]; 
(*With probability q, flip deviants*) 

(*For PSP use the below flipper instead*) 

(*Flipper [p_, average_, q_] := 

Table [RandomChoice[{l - q, q} -> ■[p[[i]], 
If [Abs [average - p[[i]]] < 0.5, 

p[[i]], Flip[p[[i]]]]}], {i, 1, u}]; 

*) 

(*Construct the initial population*) 

p[[l]] = Table [RandomChoice[-[0, 1} -> {0, 1}] 
, {i , 1 , u}] ; 

(*measure the initial average 
via chosing n random people*) 

variance [ [1] ] = Variance [RandomSample [p [ [1] ] , n]]; 

(*iterate*) 

For[i =2, i <= steps, i++, 

picker = RandomSample [p [ [i - 1]], n] ; 
(♦select n people*) ; 

average [[i - 1]] = Mean [picker] ; 
(*Measure their average value*) ; 

variance [ [i] ] = StandardDeviation [picker] ; 
(♦Measure their variance*) ; 

p[[i]] = 

Flipper [p [ [i - 1] ] , average [ [i - 1] ] , q] ; 
(*Flip everyone accordingly*) ; 
] 

grandaverage = Table [Mean [p [ [i]]] , {i, 1, steps}]; 

(*Scientific Population Simulation*) 
ensemblesize = 1000; 
timer ange = 26; 
epsilon=10"{-7} ; 
Clear [ensemble, w] ; 

ensemble = Table [0, {i, 1, ensemblesize}]; 
For[w = 1, w <= ensemblesize, w++. 

Clear [i, j, p, data, 

publishprob, error, q, x, av, t, b] ; 

steps = 90; 

timerange = 26; 

trail = 4; 

p = Table [{0, 0}, {i, 1, steps}]; 
average = Table [0, {i, 1, steps}]; 
truth = 1; 
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initiall = 0.145; 
iiiitial2 = 0.141; 
error [t_] : = 
If [t <= timerange , 
initiall (1 - t/ (timerange + 1)) .epsilon] ; 
p[[l]] = {1, 1 + initiall}; 
p[[2]] = {2, 1 + initial2}; 

q = 1; 

j = 2; 

trailingdeviation = Table [0, {i, 1, steps}]; 
trailingaverage = Table [0, {i, 1, steps}]; 
trailingdeviation [[1]] = 
StandardDeviation[{p[[l, 2]],p[[2, 2]]}]; 
trailingaverage [[1]] = (p[[l, 2]] + p[[2, 2]])/2; 

For[i = 2, i < steps, i++, 
data = REmdomReal [NormalDistribution [truth, 

error [j]]] ; 
trailingdeviation[ [j] ] = 
If [j == 1 I I trail ==1,0, 
If [j <= trail, 

StandardDeviation [Table [p[[k, 2]] , 
{k, 1, j}]], 

StandardDeviation [Table [p [ [k , 2] ] , 
{k, j - trail + 1, j}]]]] ; 
trailingaverage [ [j] ] = 
If[j <= trail. Mean [Table [p[[k, 2]], 



{k, 1, j}]], 
Mean [Table [p[[k, 2]] , 

ik, j - trail + 1, j}]]] ; 
If [Abs [trailingaverage [ [j] ] - data] <= 
q (error [j] + trailingdeviation [[j]] ) , 

p[[j]] = {j, data}. True] ; 
]; 

p = Table[p[[i, 2]] , {i, 1, j}] ; 

ensemble [ [w] ] = p 

] 

Show [ListPlot [Take [p, timerange - 1], 
PlotRange -> {0.8, 1.4}, 
Mesh -> Full, PlotStyle -> Red], 
Plot[l + (p[[l]] - 1) Exp[-(x - l)/t], 

{x, 1, steps}, 

PlotRange -> All, PlotStyle -> Red]] 
averages = 

Table [Mesoi [Table [ensemble [ [i , k] ] , 
{i, 1, ensemblesize}] ] , {k, 1, timerange}] 
errors = Table [ 

SteindardDeviation [Table [ensemble [ [i , k] ] , 

{i, 1, ensemblesize}]], {k, 1, timercoige}] 
ListPlot [averages, PlotStyle -> {Green}, 
PlotMarkers -> {Automatic, 7}] 



